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We discuss preliminary results on medium-modified fragmentation functions obtained 
in a combined NLO fit to data on semi-inclusive deep inelastic scattering off nuclei and 
hadroproduction in deuteron-gold collisions. 

1 Motivation 

In the last few years there has been a significant improvement in the perturbative QCD de- 
scription of hadroproduction processes, and more specifically, in the precise determination of 
fragmentation functions [T] . One of the most interesting features of these extractions is that 
they not only reproduce the usual data on electron-positron annihilation into hadrons, but 
they describe as well other processes like semi-inclusive deep inelastic scattering (SIDIS) and 
hadroproduction in proton-proton collisions, with remarkable precision. This demonstrates 
the ideas of factorization and universality, which are the starting points for the perturbative 
QCD description, in the kinematical domain accessed by these experiments. 

In spite of this apparently successful scheme, fragmentation functions (FFs) are far from 
giving us the whole picture of hadronization; they are an efficient way to encode the non- 
perturbative information required for computing cross sections, but they don't tell us much 
about, for instance, the space-time evolution of the hadron formation, the origin of the 
differences between the different species, and many other non perturbative key issues. 

To have an insight into these features, we need to attack the problem from a different 
angle, that could be for example, studying how hadronization proceeds in a nuclear environ- 
ment. In recent years there has been an increasing number of experiments pointing in this 
direction, exploring hadroproduction off nuclear targets with mass A, such as very precise 
semi-inclusive eA measurements by HERMES [2] , as well as experiments using dAu collisions 
at the RHIC [3l |4| . Both processes show clear signals of non trivial nuclear dependence in 
the hadronization mechanism, that needs to be understood. Additionally, this information 
serves as baseline for heavy ion research programs ongoing at RHIC and projected for LHC, 
aiming at unravelling the properties of hot and dense QCD matter [5] . 

In parallel with the availability of data, several theoretical calculations and models de- 
signed to describe the nuclear dependence have been proposed, some of them putting the 
emphasis on the interactions of hadrons and "pre-hadrons" in the nuclear medium, others 
focussing on the interactions of the "seed partons" with the medium. More recently, it has 
also been suggested that FFs could obey different energy scale dependence because of the 
nuclear medium. Most models reproduce with different success features of the data, in spite 
of very different, even orthogonal, approaches and ingredients [S]. 

In our approach, rather than proposing a particular mechanism, in a first step we try to 
isolate the medium induced modification of the FFs as precisely as possible, factoring out 
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the standard knowledge on nuclear effects in the parton densities, as seen in deep inelastic 
scattering or Drell Yan processes with nuclei [5], and the QCD dynamics inherent to hard 
processes at NLO accuracy. In other words, we fit medium-modified fragmentation functions 
(nFFs), assuming (or testing) factorization and universality in a nuclear environment. 

2 Convolutional Approach: 

Rather than fitting from scratch the nFFs, which would take as many parameters as the 
standard or vacuum FFs, plus some more to represent the nuclear size or density dependence, 
in the following we choose to relate the nFFs to the standard ones by a convolution: 

^fM(^,Qo)= f'^W,iy,A,Ql)D':{^ Ql), (1) 

J z y y 

where the weight function Wi{y, A, Qq) parameterizes nuclear effects in FFs at a given initial 
scale Qq. The scale dependence of nFFs is then determined by ordinary NLO evolution 
equations. A simple delta function (5(1 — y) as weight would imply no nuclear effects, while 
a shift in its argument, i.e. 6{1 — e — y), would represent a shift in the momentum fraction 
as suggested, for example, by some energy loss mechanism. A more general weight function 
like 

W,{y,A,Ql)=n,y'''{l-yf', (2) 

can parameterize effects not necessarily related to modifications the parton's momentum, 
such as hadron or pre-hadron attenuation or enhancement, with a great economy of parame- 
ters and, at the same time, retaining the information on the vacuum FFs. The A dependence 
of the weights can easily be included in the coefficients taking them as smooth functions 
of A. On the other hand, convolutional integrals arc the most natural language for parton 
dynamics beyond the LO and allow the application of the McUin transform technique in the 
NLO computation of the scale dependence and cross section estimates. 

The convolutional approach has been shown to be specially effective in the extraction of 
initial state nuclear effects in inclusive DIS and Drell Yan processes at NLO accuracy [6] . For 
nFFs this is also the case, and it can be shown that an extremely simple functional form for 
the weights with very few parameters can reproduce the main features of the available data. 
However, since we are interested in obtaining an accurate parameterization, we implement 
more flexible weights to account for mechanisms other than global shifts in momentum but 
still preserving factorization. Specifically, we adopt the following ansatz: 

Wgiy,A,Ql)=n,d{l-e,~y) + n'^y''-il~yf- (3) 
Wg{y, A, Ql) = ngd{l ^ ^ y) + n'^ y^^ (1 - yf^ (4) 

discriminating between quarks and gluons. In order to enhance the sensitivity of the fit 
to the gluon fragmentation, it is essential to include in the fit also hadroproduction rates 
from dAu collisions at RHIC [31 H] . The A dependence of the coefficients is implemented as 
"n-i = A„. + 7„. A^'^i . In most cases, however, A can be set to zero or to unity due to the 
vanishing of nuclear effects as A ^ 1. The assumption of a linear behavior (5 = 1) does not 
spoil the quality of the fit. 

Results of a fit with such an ansatz can be found in Fig. 1-3, for SIDIS pion multiplicities 
off different nuclei from HERMES [5] (rates are normalized to a deuterium target), neutral 
and charged pion production in dAu collisions from PHENIX and STAR [4j , respectively. 
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The dashed hnes represent 
the estimates for these cross 
sections computed by taking 
into account nuclear effects 
only for the PDFs. The soUd 
hnes refer to our results with 
nuclear modifications also in 
the FFs. We obtain an over- 
all = 350.45 for 368 data 
points and 17 parameters, re- 
sulting in x^/d.o.f = 0.997. 
Focusing on SIDIS data in 
Fig. 1, the quality of the fit 
is quite impressive, and even 
the XBj (or v) dependence 
is reproduced although the 
weights in Eqs. 3 and 4 and 
the nFFs only depend on z as 
required in a factorizcd ap- 
proach. In addition, there 
seems to be no conflict with 
the standard (vacuum) 
dependence assumed in our 
fit. RHIC data are also re- 
produced well as shown in 
Figs 2 and 3. 
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Figure 1: Pion SIDIS multiphcitics from HERMES 
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Figure 2: tt*^ production in dAu collisions 



The logarithmic fall off of the cross sec- 
tion over four orders of magnitude hinders 
the comparison between data and the the- 
oretical estimates with and without nFFs. 
The differences are sizable, as shown in the 
lower panel of Fig. 2: filled and open squares 
are the results for "(data-thcory)/thcory" 
for fits with and without modified FFs, re- 
spectively. PHENIX data indicate a signif- 
icant suppression at larger pT , while STAR 
data only hint at such a behavior towards 
the end of the more limited px range. Even 
though the fit follows the data down to 
very low px, this needs to be taken with 
extreme caution since possible soft contri- 
butions may spoil factorization, and stan- 
dard PDFs and FFs do not give a satisfac- 
tory description of pp data at NLO below 
PT =i 2 GeV. 
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In Table 1 we present 
the values obtained for the 
17 free parameters of the 
fit. Table 2 shows the values 
the coefficients in the weight 
functions in Eqs. 3 and 4 take 
for different nuclei using the 
parameters in Tabic 1. As 
it can be noticed, the A de- 
pendence of the coefficients is 
rather smooth and typically 
linear for A not too small. 
Quark and gluon coefficients 
are found to be rather differ- 
ent. In the case of quarks 
(and antiquarks) the main 
modification due to the nu- 
clear medium is related to 
the first term in Eq. 3, asso- 
ciated to a shift eq in the mo- 
mentum fraction. This shift 
is typically small, but introduces a non trivial z dependence, and increases linearly with A. 
The normalization of this term, in turn, decreases with A leading to an overall suppression. 
More than 80% of the nuclear modification comes from this term, while the second term 
represent a smaller effect. For giuons, both the first and the second term compete resulting 
in a much more pronounced effect. 
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Table 1: Parameters describing 


the medium-modified FFs at NLO 
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Table 2: Coefficients of the weight functions in Eqs. 3 and 4 for different nuclei 



It is also enlightening to analyze the ratios of the nuclear and vacuum fragmentation 
functions as a function of z for a given Q^, as shown in the left panels of Fig. 4. We also 
show the nFFs themselves in the right panels, and always at = lOGeV^. 
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For quarks and anti- 
quarks the main effect seems 
to be a reduction in the 
fragmentation probabihty in- 
creasing with z and with nu- 
clear mass A. There is a sig- 
nificant drop at smah z (z < 
0.1), however, neither the nu- 
clear nor the vacuum frag- 
mentation functions are well 
constrained by the data, and 
neglected final state hadron 
mass effects can be signifi- 
cant. Gluon fragmentation 
shows a rather different pat- 
tern of medium modification, 
with a much more noticeable 
drop at large z, combined 
with a sizable enhancement 
at intermediate z. Notice 
that the low z behavior of 
the gluon nFFs is crucial for 
reproducing the XBj depen- 
dence shown by the SIDIS 
rates and therefore it is rel- 
atively well constrained by 
these data. On the contrary, 
the large z behavior is re- 
lated to the suppression of the RHIC cross sections at low that need better understanding 
from the theoretical side. 




Figure 4: nFFs and rates at = 10 GeV^ 



3 Conclusions 

We have explored the feasibility of reproducing the main features of data sensitive to medium 
induced modifications of fragmentation functions in a factorized QCD framework at NLO 
accuracy based on a convolutional approach. The nice agreement found encourages us to 
further pursue this approach and extend it to hadrons species other than pions. 
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